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REVISION 1 CHANGES 


This document represents a substantial revision of the previous one dated 
October 1971. The major modification involves the use of the Unified Powered 
Plight Guidance Routine during the deorbit targeting to accurately estimate the 
effects of the finite thrust maneuvers for the low thrust/weight ratios associated 
with the Shuttle project. 
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FOREWORD 


This document is one of a series of candidates for inclusion in a future 
revision of JSC -04217 "Space Shuttle Guidance, Navigation and Control Design 
Equations". The enclosed has been prepared under NAS9- 10268, Task No. 15- A, 
"GN&.C Flight Equation Specification Support", and applies to functions 1, 2, 4 
and 5 of the Entry Guidance Module (OG5) as defined in JSC -03690, Rev. D , 
"Space Shuttle Orbiter Guidance, Navigation, and Control Software Functional 
Requirements", dated January 1973. 
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1 . 


INTRODUCTION 


The large entry crossrange capability of the Shuttle permits deorbit to a 
specified landing site to be accomplished with a single maneuver. Since the re- 
quired velocity change is smallest when no plane change is made, the equations 
presented here are designed to target the Unified Powered Flight Guidance Routine 
(Ref. 1) for an in-plane maneuver. The ignition time for this maneuver is selected 
to satisfy entry interface and landing site constraints with minimum fuel expendi- 
ture. 

If the Shuttle had no crossrange capability, an in- plane deorbit to a speci- 
fied landing site would occur only when the landing site, which rotates with the 
earth, intersects the orbital plane of the vehicle. Assuming the landing site lati- 
tude is less than the orbital inclination and neglecting the effects of precession, 
the landing site will intersect the orbital plane twice every twenty-four hours. The 
time difference between these two intersections is, in general, not twelve hours. 

In the case when the landing site latitude is equal to the orbital inclination there 
will be only one intersection every twenty- four hours. 

Since the Shuttle has a large crossrange capability, deorbit can occur when- 
ever the angle between the landing site vector and the orbital plane is less than 
approximately 20 deg. In general, there will be two sets of opportunities every 
twenty-four hours. Within each set, there may be several deorbit opportunities 
occurring on consecutive orbits with varying crossrange requirements. When the 
latitude of the landing site approaches the inclination of the orbit, these two sets 
merge to become one. If the landing site latitude is greater than the orbital in- 
clination, the landing site may still fall within the crossrange capability of the 
vehicle. Based on these considerations, this routine has been designed to step 
through successive solutions, allowing the crew to select a particular deorbit 
opportunity based upon entry crossrange, time -to -ignition, required velocity change, 
landing site lighting conditions, urgency of the return, etc. 

The deorbit targeting program is divided into two sections. The first sec- 
tion establishes an initial guess of the landing time based on updating the vehicle in its 
orbit until its position coincides with the landing site position projected into the plane 
of the orbit. If the crossrange at this point is not acceptable, the vehicle is updated 
through successive orbital periods until an acceptable crossrange distance is at- 
tained. The second section utilizes the Unified Powered Flight Guidance (UPFG) 
routine - which is based on velocity-to-be-gained and Linear Tangent Guidance 
concepts — in order to approximate the effects of the finite deorbit maneuver. This 
is necessary as the anticipated thrust/ weight ratio for the Shuttle during deorbit 
using a single OHMS engine is 0.025, resulting in long burn arcs. The UPFG rou- 
tine is called in an iterative loop \vhich searches for the minimum fuel deorbit trajec- 
tory. 


l-l 



The desired entry interface conditions (at 400, 000 ft) are assumed specified 
by the relationship 

v r = c l + c 2 v h 

where v and are the radial and horizontal components of velocity at the entry 
interface point and and Cg are referred to as the entry interface parameters. 
This condition is approximately equivalent to a linear relationship between the entry 
interface angle and velocity. The flight time and range between entry interface 
and the landing site are assumed to be functionalized in terms of the entry interface 
flight path angle. 
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2 . 


FUNCTIONAL FLOW DIAGRAM 


A functional flow diagram presenting the basic approach to the deorbit 
targeting problem is presented in Figure 1. In addition to the vehicle's state 
vector, the major inputs are landing site location (latitude and longitude) , entry 
interface parameters, earliest desired time of landing and maximum cross range. 

The first section in the deorbit targeting program is to select the earliest 
landing time which results in an acceptable crossrange. This step is computed 
using the constraint that the vehicle must arrive above the in- plane projection of the 
landing site at approximately the same time as was assumed in determining the 
location of the landing site (which rotates with the earth). In this first step, an 
estimate of the difference in the time flight between (1) the actual time spent 
from deorbit to landing and (2) the orbital time associated with the corresponding 
central angle is utilized. This estimate At p can be represented by a constant 
as the following steps in the program will compensate for any errors due to this 
assumption. 

To accomplish the first step, the state vector is precision updated to the 
time tp L g (initially set' approximately equal to the desired earliest time of landing). 
After obtaining the landing site's inertial position based on the time t pj g minus 
At^, , the time Atj p associated with a conic update of the state vector to a point 
colinear with the projection of the landing site into the orbital plane is computed. 

A new estimate of the landing time is now obtained by adding the time At ip to the 
previous landing time. After precision updating the trajectory to the time t pLg 
plus At ip , another calculation of Atj p is made. Using this At jp in a new 
calculation of the landing time, the above constraint is assumed satisfied. This 
procedure is illustrated in the following sketch. 




If the crossrange d CR computed by the above procedure is not acceptable, the 
landing time is increased by an orbital period and the above procedure is repeated - 
with the exception that the second conic update is bypassed - to establish a new 
crossrange. This process is repeated until an acceptable crossrange is attained. 

The second section starts with a computation of an estimated deorbit time 
and an entry interface position. Following a conic update of the state vector to the 
deorbit time, an approximate computation of the optimum-inertially defined-turning 
rate for the thrust direction during the burn is made. The Unified Powered Flight 
Guidance Routine is next called to compute the variable • This variable 

represents the approximate direction of the thrust vector at the midpoint of the ma- 
neuver, measured with respect to the local horizontal direction. It has been shown 
(Refs. 2-6) that a burn arc which utilizes a thrust direction which rotates at a con- 
stant inertial rate and which results in a value of zero represents an ap- 

proximate fuel optimum maneuver. The fuel optimum deorbit trajectory can there- 
fore be computed by iteratively varying the time of the deorbit maneuver until 
o'mid i- s driven to zero. This is accomplished using the Newton-Raphson iterative 
technique incorporated in the ITER subroutine (Ref. 7). The iteration is accomp- 
lished in two phases. In the first phase, the entry interface position and landing 
time are held fixed while the orbit is conically updated to the ignition point on each 
iteration cycle. In the second phase - when convergence is near - the entry interface 
point is computed based on the latest estimate of the landing time and the orbit is 
precision updated to the ignition point on each iteration cycle. During this phase, 
the UPFG routine computes a deorbit maneuver based on a precision update of the 
trajectory from the end of the burn to the entry interface point. 


Ref. 1, containing the description of the UPFG routine, does not presently 
allow for this precision updating mode. It will, however, be included in 
a future update of this reference. 
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DEORBIT TARGETING 



Figure la. Functional Flow Diagram 
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Figure lb. Functional Flow Diagram 
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3. INPUT AND OUTPUT VARIABLES 

The inputs to the deorbit targeting program are; 

£ O’ — 0 State vector 

Time associated with r Q J v Q 
tgTL Earliest desired time of landing 

Latest desired time of landing 
m Estimated vehicle mass 

s OMS 1 Number of orbital maneuvering system engines 

<P LS Landing site latitude 

X LS Landing site longitude 

^ACR Maximum acceptable crossrange distance 

The output variables are: 
t. Ignition time 

tpnj Time of arrival at the offset target vector 

Ay Maneuver Ay 

r Qp Offset target vector 

n rev Integral number of revolutions 

Jijvj Unit normal to transfer plane in direction of angular momentum vector 

s goln Lambert solution switch 
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4. 


DESCRIPTION OF EQUATIONS 


The majority of the deorbit targeting program consists of calls to other 
routines and simple logical operations, interspersed with some elementary equa- 
tions. Those operations which are not self evident in terms of the functional flow 
description and the detailed flow diagram are discussed below. 


4 . 1 Computation of 8 ^ p 

The angle between the position vector r ^ specified at the projected 
landing site time t p ^g and the projected landing site vector _i.^g p is defined 
as 8 ip • Whenever the crossrange computed in the first section of the deorbit 
targeting program exceeds the maximum acceptable crossrange, the vehicle's 
state vector is updated through one orbital period in an attempt to establish an ac- 
ceptable crossrange. During this interval of time, the inertial position of the 
landing site changes (due to the rotation of the earth), resulting in an in-plane 
projection of the landing site which either proceeds or lags the updated state vec- 
tor. An upper limit on this difference (0.j. p ) can be established as +30° . When 
the rotational direction of the orbit is opposite to the rotational direction of the 
earth, this difference will be a negative quantity. Hence, 0j. p can be computed 
as an angle between the limits of -30° and 330° . 


In order to make the first computation in Section 1 compatible with sub- 
sequent computations, the state vector is initially advanced by one -twelfth of an 
orbital period (approximately equivalent to 30° ) beyond the earliest desired time 
of landing t pTL • 


4. 2 Computation of a ) p 

Based on the primer vector theory, Refs. 2, S and 9, the thrust direction 
cx during an optimum burn must be specified by 


tan ex - -sin 0 / (c^ + 2 cos 9 ) {4, 1) 

0 is the angle from the burn centroid to the point where a is desired. This for- 
mulation assumes that the burn occurs under near circular orbital conditions with 
an inertially defined- constant turning rate maneuver. The constant can be 
obtained from the entry interface constraint that the primer vector must be per- 
pendicular to the entry interface constraint line, whose slope is given by (see equa- 
tion in Section 1) 

d v r / d v h = c 2 . (4. 2) 
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The desired direction for the primer vector is 

flp = - [tan Ncg) + 7r / 2 ] (4. 3) 

As the thrrst direction coincides with the direction of the primer vector, Eqs. 

(4. 1) and (4. 2) can be combined with the result 

c 3 = “ 2 cos ® ME1 " c 2 sin ® ME I ^ 4 ‘ 4 ^ 

where 0 mei an S^ e between the burn centroid and the entry interface point. 

Equation (4, 1) now defines the optimum thrust angle as a function of 0, 
®MEI anc * c 2 * Using a central angle © M which extends from the centroid to 
the beginning of the maneuver { 0 M is negative), the thrust direction at 

the start of the maneuver can be computed. As the thrust direction at the burn 
centroid should be zero {see Eq, (4. 1)), the inertial ly defined turning rate Wp 
for the thrust direction is obtained by noting that the thrust vector must be rotated 
through both the angle and 0 ^ 

(Up = + (4.5) 

where tg is the estimated burn time. 

The above computations are performed once for each iterative cycle in the 
second section of the program. Initially, (Op will be only approximately correct 
due to an inaccurate estimate of Av , As the iteration converges to an acceptable 
value of Q'n^jjp > the computation of (Op should become stable and sufficiently 
accurate (the computation of is not very sensitive to changes in Wp so the 

approximation involved in the computation of <Up are inconsequential). 
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5 . 


DETAILED FLOW DIAGRAM 


This section contains detailed flow diagrams of the deorbit targeting rou- 
tine. 

Each input and output variable in the routine and subroutine call statements 
can be followed by a symbol in brackets. This symbol identifies the notation for 
the corresponding variable in the detailed description and flow diagrams of the 
called routine. When identical notation is used, the bracketed symbol is omitted. 
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Figure 2b. Detailed Flow Diagram 
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Call Iteration Routine (Ref. 7) 
Input: c, [©]» "tp [ x ]> 

ff MID0 f e 0 J^DO^O J' 4 
Output: C, tjp [x], f^Q ]» 


*00 fx 0 ^ s fail 



Figure 2d. Detailed Flow Diagram 
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Figure 2e. Detailed Flow Diagram 
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6 . 


SUPPLEMENTARY INFORMATION 


The deorbit targeting technique presented here requires long term state 
vector extrapolation to select a deorbit opportunity several revolutions later. The 
accuracy of this extrapolation is dependent upon accurate Knowledge of the current 
state vector. Any errors due to imperfect navigation will tend to be magnified by 
the long term extrapolation. Fortunately, the out- of- plane component of error 
tends to oscillate, and can be expected to remain below about 500 ft. Consequently, 
the prediction of crossrange for a deorbit opportunity several revolutions later 
is not significantly affected by the expected state vector error. Mission planning 
functions, which are sensitive to crossrange requirements, can be carried out 
without being significantly affected by the navigation error. 

The in-plane component of state vector error does increase during long 
term extrapolation, at a rate of about 1000 ft per revolution. Therefore the tar- 
geting should be repeated during the revolution prior to the deorbit maneuver , 
taking advantage of more precise knowledge of the vehicle state vector. 

The effect of different azimuths on the entry phase has not been considered 
in this design. Since the crossrange capability of the vehicle is dependent on 
azimuth, and thus may be larger in one direction than the other, this should be 
considered in the acceptable cross range criterion. 

The equations for the determination of the time and range from entry inter- 
face to landing are not yet available. Additional analysis in the area of entry 
guidance should allow the functionalizations of these parameters in terms of the 
entry interface flight path angle. 

The de orbit targeting equations presented here are designed to target a 
single minimum fuel maneuver which satisfies entry interface and landing site con- 
straints. Since deorbit is the final major maneuver, it may be logical to use all 
the remaining fuel to either effect a faster return, to reduce the entry cross range 
requirement or to decrease the deterioration of the heat shield (by lowering the entry 
interface velocity). Various methods for accomplishing these goals may be used , 
including thrusting out-of-plane, using multiple maneuvers to place the vehicle in 
a phasing orbit prior to the deorbit maneuver and minimizing the entry interface 
velocity in place of minimizing the fuel. For the most part, these modifications 
to the deorbit targeting would be superimposed on the logical structure herein 
contained. 

Various empirical studies were undertaken to verify the new deorbit targeting 
logic. These studies are covered in Refs. 2-6 and include data for circular orbits 
between altitudes of 70 n. mi. and 500 n. mi. and for a 100 x 270 n. mi. elliptical 
orbit. 


6-1 



REFERENCES 


1. Brown, D. , Higgins, J. , Brand, T, , "Unified Orbiter Powered Flight 
Guidance 11 , Space Shuttle GN&C Equation Document No. 24, Draper 
Laboratory, June 1973. 

2. Tempelman, W. , "A Survey of Fuel Optimum Impulsive and Near Fuel 
Optimum Finite Thrust Maneuvers for Deorbits from Circular Orbits 
and a 100 x 270 n. mi. Elliptical Orbit", Draper Laboratory, 23A STS 
Memo No. 28-73, 29 March 1973. 

3. Tempelman, W. , "Deorbit Trajectories From a 100 x 270 n. mi. Ellip- 
tical Orbit Based on Fixed Deorbit Position”, Draper Laboratory, 23A 
STS Memo No, 31-73, 6 April 1973. 

4. Tempelman, W, , "A Survey of Deorbit Trajectories Using a Target 
Line Constraint", Draper Laboratory, 23A STS Memo No. 34-73, 

23 April 1973. 

5. Tempelman, W. , "Deorbit Targeting Using the Unified Powered Flight 
Guidance Program”, Draper Laboratory, 23A STS Memo No. 37-7 3, 

27 April 1973. 

6. Tempelman, W. , "Comparison Between Predicted Deorbit Targeted 
Maneuvers and the Actual Maneuvers”, Draper Laboratory, 23A STS 
Memo No. 40-73, 8 May 1973. 

7. Tempelman, W. , "Rendezvous Targeting”, Space Shuttle GN&C Equa- 
tion Document No. 7, Rev. 2, Draper Laboratory, February 1972. 

8. Brown, D. , "Thrust Direction Turning Rate for Large Central Angle 
Shuttle Deorbit Burns", Draper Laboratory, 23 A STS Memo No. 55-73, 
15 June 1973, (to be published) 

9. Lawden, D. , "Optimal Trajectories for Space Navigation”, Butterworths, 
London, 1963, p. 96. 

10, Shepperd, S. , Robertson, W. , "Conic State Extrapolation”, Space 
Shuttle GN&C Equation Document No. 25, Draper Laboratory, June 
1973, (to be published). 

11, Robertson, W. , "Precision State and Filter Weighting Matrix Extrapola- 
tion”, Space Shuttle GN&C Equation Document No. 4, Rev. 3, Draper 
Laboratory, October 1972. 


R-l 



